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Abstract. We investigate the effect of pair creation on a shock struc- 
ture. Actually, particles accelerated by a shock can be sufficiently ener- 
getic to boost, via Inverse Compton (IC) process for example, surrounding 
soft photons above the rest mass electron energy and thus to trigger the 
pair creation process. The increase of the associated pair pressure is thus 
able to disrupt the plasma flow and possibly, for too high pressure, to 
smooth it completely. Reversely, significant changes of the flow velocity 
profile may modify the distribution function of the accelerated particles, 
modifying consequently the pair creation rate. Stationary states are then 
obtained by solving self-consistently for the particle distribution function 
and the flow velocity profile. We discuss our results and the application 
of these processes to the high energy emission and variability of compact 
objects. 



1. Introduction 

The high energy emission observed in compact objects like AGNs or X-ray bi- 
naries requires the existence of high energy particles. Shocks being particularly 
attractive particle acceleration sites, they are generally believed to occur in the 
central region of these objects. In these cases however, the presence of important 
radiation fields necessitate to take into account particles-photons interactions in 
the shock region. 

If some works have already studied particles acceleration at shocks including 
radiative cooling processes (Webb et al., 1984; Drury et al., 1999) we have in- 
vestigated the effect of pair creation on the shock structure (Petrucci et al., 
hereafter POO). Particles accelerated by the shock can effectively be sufficiently 
energetic to boost, via Inverse Compton (IC) process for example, surrounding 
soft photons above the rest mass electron energy and thus enable to generate 
pairs. Consequently, the increase of the pair pressure may be able to modify the 
plasma flow and eventually, for too high pressure, to smooth it completely. 
We describe here the geometry of the (toy) model we use (cf. Fig 1). We present 
the basic equations and the main results of this model. We then briefly discuss 
its application to the high energy emission and variability of compact objects. 
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2. The Toy Model 

The schematic view of our toy model is plotted in Fig. 1. We suppose the exis- 
tence of a thermal supersonic plasma undergoing an adiabatic non-relativistic 
shock. It is supposed to be embedded in an isotropic external soft pho- 
ton field. We assume the presence of a magnetic field B D , parallel to the shock 
normal, frozen in the plasma and slightly perturbed by Alfven waves. We sup- 
pose the magnetic energy to be in equipartition with the particles kinetic energy 
so that, in each part of the shock, the Alfven wave speed is equal to the flow 
speed. Particles are then scattered by Alfven waves through pitch angle scat- 
terings, gaining energy through the first order Fermi process when crossing 
back and forth the shock discontinuity. We also suppose the magnetic perturba- 
tions to have sufficiently small amplitudes so that we can treat the problem in 
quasilinear theory using the Fokker-Planck formalism. Finally, we suppose 
the plasma pressure to be dominated by the pressure of the relativis- 
tic particles. In this case, the flow velocity profile makes already a smooth 
transition between the up and downstream region but acceleration still occurs. 



2.1. The Geometry 



The shock location: We assume a ID geometry. It means that the different 
parameters characterizing the flow are supposed to be homogeneous in each 
section perpendicular to the x axis. This is a relatively good approximation in 
the central parts of the shock where the border effects are negligible. We suppose 
the shock to be located in x = (cf. Fig. 1). During numerical integrations, this 
will be ensuring by imposing that the flow velocity u(x) possesses an inflection 
point in x = that is: 

=0 (1) 

=0 



dx 2 



The acceleration region: A particle of Lorentz factor 7 will interact with the 
shock if it is located within about one diffusion length £.0(7) = D xx /u (where 
D xx is the spatial diffusion coefficient) from the shock. On the other hand it 
will be cooled (by Inverse Compton process) on a cooling length scale L coo i(~f). 
Ld and T coo i are increasing and decreasing functions of 7 respectively. There 
thus exists a Lorentz factor j c for which: 

L D {lc) = L coo i(j c ) = L 

We will define the acceleration region as the physical space —L<x<L (hatched 
region in cf. Fig. 1). Consequently, in this region the coolings are negligible 
in comparison to heatings for particles with 7 < 7 C since for such particles 
L D (j) < L cooi{l)- 

In the following we will suppose the spatial diffusion coefficient D xx to be 
independent of the energy of the particles and it will be simply written D. 
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Figure 1. The geometry of our (toy) model. The shock discontinuity 
is represented by the bold line. The hatched region is the acceleration 
region and the pair creation one is filled in pink. Particles are rep- 
resented by straight arrows and photons by warped ones (in blue for 
surrounding soft, i.e. UV, photons, and in red, high energy, i.e. X-ray, 
photons). Scales are not respected. 



The pair creation region: Particles accelerated in the shock will produce high 
energy photons by scattering, via Inverse Compton process (IC), the external 
soft photons. These high energy photons will interact with themselves to produce 
pairs on a length scale i? 77 . Typically we have: 

R — 

1 + T 77 

where R is the typical size of the shock (cf. Fig. 1). This expression take into 
account the photon-photon depletion for r 77 3> 1 and the geometrical dilution 
of the X-ray photon density, far from the shock, which upperlimits the pair cre- 
ation region size to roughly the shock size R. We will define the pair creation 
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region by — i? 77 < x < +iiy 7 . 

We assume the pair pressure creation rate to be constant in the pair creation 
region (that is for — i? 77 < x < -R 77 ) and null outside. This assumption may 
be supported by the fact that magnetic turbulence, supposed to be present to 
scatter particles, is a useful process to isotropize and homogenize the plasma in 
the vicinity of the shock in very small time scale 

The "effective" compression ratio: The definition of the compression ratio like 
the ratio between the far upstream and downstream flow velocity is rather un- 
satisfactory in our case since it will not give a real estimate of the velocity 
change experienced by a relativistic particle in the vicinity of the shock, where 
the particle distribution is principally built. We will thus define an "effective" 
compression ratio, noted simply r, as the ratio between the upstream velocity in 
— 10L and the downstream velocity in +10L. We have checked that in the case 
of a strong shock without pairs, this definition still gives a compression ratio 
very near the expected value (in a plasma dominated by relativistic particles) 
of 7. 

2.2. Important energetic al thresholds 

The acceleration threshold: Only particles having a Larmor radius compa- 
rable to the wavelength of the Alfven spectrum will undergo scatterings ( Jokipii, 
1976; Lacombe, 1977) and thus will go back and through many times across the 
shock. In — p + plasma as those we deal with, the non-relativistic protons 
limit the Alfven waves spectrum to wavelength greater than 2ttVa/^cp ( ^cp is 
the cyclotron pulsation of the protons in a magnetic field B, and Va is the Alfven 
velocity). Consequently there exists a lower Lorentz factor for a relativistic lep- 
ton to be accelerated in a shock: 

_ rrip V A 

7 ^ Trrtin — 

m e C 

In the more general case "y m in is at least of the order of a few. We will admit that 
pre-accelerator processes exist (like magnetic reconnection, whistler) to bring 
particles above j m in- Consequently, since particles annihilate preferentially for 
7 ~ 1 (Coppi Sz Blandford, 1990), we will neglect the annihilation process 
in the shock region. 

The pair creation threshold: We will suppose for simplicity that the external 
soft photon field is mono-energetic with a mean photon energy e s (in m e c 2 unit). 
Thus, on average, a soft photon scattered by a lepton of Lorentz factor 7 will 

be boosted, via IC, to an energy e ~ — 7 2 e s (Rybicki & Lightman, 1979). The 

3 

high energy photon produced will be able to generate a pair electron-positron if 
at least: 

/3m e c 2 \ 1/2 
e > 2m e c 2 i.e. 7 > I ^ J = j th 

We assume j t h > Imin so that particles have to be accelerated in the shock to 
initiate the pair creation process. 
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2.3. The particle distribution 

The spectral index: It can be shown that the solution of the evolution equation 
of the particle distribution function including the pair creation process (cf. Eq. 
(5)), still has a energy power law dependence, as it is effectively the case without 
pairs (POO). The spectral index s (here s is the spectral index of the spatially 
integrated distribution function of the particles n{^f) oc 7 _s ) keeps also the same 
expression in function of the compression ratio i.e.: 

a(r) = (r + 2)/(r-l). (2) 

For plasma dominated by relativistic pressure, the compression ratio is neces- 
sarily smaller than 7 so that s(r) is larger than 1.5. 

The high energy cut-off: Since we take into accounts the cooling, a high energy 
cut-off in the particle distribution must necessarily appear at a Lorentz factor 
7 C where heating and cooling balance (Webb et al., 1984). Since we suppose 
that particles cool via inverse Compton process (assuming that the external soft 
photon density is homogeneous in the shock region) and that they are accelerated 
by the first order Fermi process, the maximum Lorentz factor 7 C achievable by 
the acceleration process may be written as follows (POO): 



1 "I r - 1 . . 



where l s is the soft compactness (l s = ?£p—l- ; me being the electron mass and 
L so f t the external soft radiation luminosity) and u\ the upstream flow velocity. 

From these different remarks, we will assume that the particle distribution 
has a cut-off power law shape, i.e.: 

n(j) oc 7 _s exp ( -— ) (4) 



1, 



3. The basic kinetic equations 

As previously said, we suppose the existence of a magnetic field B , perpendic- 
ular to the shock front and slightly perturbed by Alfven waves. The particles 
are thus scattered by these waves trough pitch angle scattering (Jokipii, 1976; 
Lacombe, 1977) and can cross the shock front several times before escaping un- 
less they are rapidly cooled by radiative processes. During these scatterings, 
the particle gain energy trough the well known first order Fermi process. We 
also suppose the magnetic perturbations to have sufficiently small amplitudes 
so that we can treat the problem in quasilinear theory using the Fokker-Planck 
formalism. Besides we assume that, in each part of the shock front, the scatter- 
ing is sufficient for the particle distribution to be nearly isotropic. With these 
different assumptions, and when first order Fermi process just as radiative losses 
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and pair creation/annihilation are taken into account, the particles distribution 
function f(p,x) must verify the following equation : 

df_ df_ ldu df_ 1 dbp A f d_ D dl 
dt dx 3 dx^ 'dp p 2 dp dx dx 

+C±(/) + D±(/). (5) 

The three first terms of the right member correspond to the first order process, 
the radiative losses (b being > 0, the detailed expression of b in the case of 
Inverse Compton cooling can be found in POO) and the spatial diffusion respec- 
tively, C±(/) is the pair creation rate and D±(/) the annihilation one. This 
equation is essentially the equation of the cosmic-ray transport originally given 
by Parker (1965), Skirling (1975 and references therein) except for the addition 
of the radiative losses and the pair processes. 

Since we suppose that the shocked plasma pressure is dominated by the 
pressure of the relativistic particles, we can deduce, from Eq. 5, the hydrody- 
namic equation linking pairs (through their pressure) and the flow velocity u(x), 
that is (in stationary state): 

dP rp] 4 du d n d „ • • 
U -ftT + 3^ = d-x D -d~x P ^ + Q «* + P± (6) 

where Q re l i s the pressure loss rate due to radiative losses and P± the pair pres- 
sure creation rate. 



A way to solve Eq. (6) is to integrate it between — oo and +L. In this case, 
we can neglect the cooling during the integration. Indeed, even if for x < —L 
particles do not interact with the shock, we have assumed that some processes 
apply and balance coolings so that particles are injected above the resonance 
threshold 7 m j n . On the other hand, in the shock region (that is — L < x < L), 
we have seen in section 2.1 that the coolings are negligible, in comparison to 
heatings, for particles with a Lorentz factor smaller than j c (the case of the 
majority of the particles). Besides, we can also neglect the annihilation since 
it occurs mainly at low energy, i.e. for particles with Lorentz factor 7 ~ 1 (cf. 
Coppi & Blandford, 1990), whereas we assume 7 > ^y m in- In these conditions, 
and with the assumption that the pair creation is homogeneous in the region 
— .Ry-y < x < +i? 77 (meaning P± roughly constant) and null outside (meaning 
P±=0), the integration of Eq. (6) between —00 and L, combined with the 
momentum conservation equation, gives (in reduced units): 



du 
dx 



I (l - a) G- a ) +n ( l+ t) (7) 



where the reduced variables are defined as followed: 



u = — 8 

Ul 



j and i? 77 = ^p- (9) 
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n = 



P±R 



77 



(10) 



II is the ratio of the pair luminosity (the pair power density integrated on the 
"ID" volume i? 77 ) to the kinetic energy flux of the flow. This parameter will 
play an important role in the evolution of the shock profile as we will see in the 
following. We can anticipate that a large value of II will be certainly unfavor- 
able to the formation of the shock. It means that all the kinetic energy of the 
upstream flow will be dissipated in pair creation processes. 



4. Shock disappearance 



_ o. 



>> 0.6 




shock 
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Pair luminosity n 



0.2 



Figure 2. Plot of the 3 solutions of Eq. (11) versus II for i? 77 = 50. 
The flow velocity tt sc h oc k at the shock location, i.e. in x = is plotted 
in dashed line. There is a maximal value of II, H m ax, above which 
only one real solution exists. Consequently, for IT > Il max the shock 
disappears. 

The shock still exists as long as Eq. (1) is verified. Combining with Eq. 
(7), Eq. (1) becomes: 

= P(u) = A9u 3 - 84u 2 + (39 + 4211)2 

x=0 

-2411 + 18^-4 = (11) 

Since P(u) is a third degree polynomial, it possesses in general 3 real solutions 
which depends obviously on II and i? 77 . They have been plotted in Fig. 2 in 

function of IT for i? 77 =50. Of course, for IT going to zero, the three branches of 
solution converge respectively to the well-known results «=l/7, 4/7 and 1 cor- 
responding to the values without pair creation. By continuity, the flow velocity 
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at the shock location will follow the second branch noted u s hock ° n the figure 
(plotted in dashed line). It appears that for a given value of i? 77 it exists 
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Figure 3. 
by ~ 0.2. 



Curve H m ax vs. Fly-/. We see that H ma x is upperlimited 



a maximal value II max of II above which there is only one real solution 
which still verifies Eq. (11). The unphysical discontinuity of the flow 
velocity at the shock location for II = Ii ma x means simply that the 
shock can not exist anymore. 

The transition between 3 to 1 real solution of Eq. (11) happens when the 
following conditions are satisfied: 



P(u) 
P'(u) 








(12) 
(13) 



The resolution of this system of equations gives thus a relation between U r 
and i? 77 : 

2/3 



1 

14 ~ 14 



63- 



IL 



R 



(14) 



77 



We have reported the corresponding function U max (R yi ) in Fig. 3. We see that 
n max is a increasing function of i? 77 meaning that the larger the pair creation 
region the larger the pair power we need to kill the shock. However, we see from 
Eq. (14) that H ma x is necessarily smaller than 3/14~ 0.20 meaning that at 
most 20% of the kinetic energy flux of the upstream flow transformed 
in pairs is sufficient to suppress the shock discontinuity. 

We have reported in Fig. 4, different shape of the polynomial P(u) and the 
corresponding flow velocity profiles obtained numerically by solving Eq. (7) for 
different values of II (fixing i? 77 to 10). We have also plotted in Fig. 5, the 
variation of the compression ratio r (as defined in section 2.1) in function of the 
pair creation rate II. As expected, r converge to ~ 1 when II increases meaning 
that the acceleration becomes less and less efficient. 
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Figure 4. disappearance of a shock due to pair creation. We have 
taken i? 77 =10 meaning that II must be smaller than ~ 0.15 (cf. Fig. 
3). 

5. Stationary states 

Up to now, we have supposed that the deformation of the flow profile, due to 
the pair pressure, does not modify the pair creation rate itself. This is a crude 
assumption since a change in the velocity profile of the flow will modify the dis- 
tribution function of the accelerated particles in such a way that the number of 
particles enable to trigger the pair creation process (i.e. particles with 7 > 7^) 
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Figure 5. Variation of the compression ratio r in function of the pair 
luminosity II for i? 77 = 1000. 



may change, modifying consequently the pair creation rate. 
From the previous section we have seen that the compression ratio r decreases 
when the pair creation rate raises (cf. Fig. 5) but reversely if r decreases, the 
spectral index will increase and the number of high energy particles will de- 
crease, the final effect being a decrease of the pair creation rate. Consequently, 
we may expect the system to reach, in some conditions, stationary states where 
hydrodynamics and pair creation effects balance. 



5.1. The parameter space 

We have studied the stationary states of our toy model by solving self-consistently 
for the particle distribution function 71(7) and the flow velocity profile u(x), the 
two function 71(7) and u(x) being linked through the pair pressure. In station- 
nary states, the system depends on six different parameters: "f m in, the minimal 
Lorentz factor for the particles to be accelerated in the shock (cf. section 2.2), 
e s , the external soft photon energy (in unit of m e c 2 ), u±, the upstream flow ve- 
locity, R the transverse size of the shock (in unit of the diffusion length L), l s , 
the compactness of the external soft photon field (cf. section 2.3) and the 

kinetic compactness defined as l^in = — 1 — ^ r . The larger l s and the larger 

AirRm e cy 

the cooling of the particles, whereas the larger lki n and the larger the kinetic 
energy of the upstream flow. 

The system is solvable in only some part of the parameter space. It always 
possesses two solutions: a "pair dominated" (large pair luminosity II, small com- 
pression ratio r) and a "pair free" one (small II, large r). The latter connects 
to the trivial solution of the problem i.e. II = 0. 
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5.2. High energy spectra 

For a given set of parameters, the system may reach stationary states charac- 
terized by a compression ratio r. The spectral index and cut-off of the particle 
distribution function (Eq. (4)) are then given by Eq. (2) and (3) respectively. 
Since we have assumed that these particles are cooled by Inverse Compton effect 
and that their distribution function follows Eq. (4), the emitted energy spec- 

r ( E \ * 

trum is characterized by a cut-off power law shape Fe oc E a exp — y— J 
where a and E c are simple functions of s and 7 C (Rybicki & Lightman, 1979): 

a = S - (15) 
E c ~ ^l 2 c e s m e c 2 (16) 




Figure 6. Contour plots of the spectral index a (solid lines) and the 
high energy cut-off E c in keV (dashed lines) of the emitted spectrum 
in the (l s , Itin) space. The soft photon energy is equal to 10 eV and 
600 eV in the left and right plot respectively. The other parameters 
have been fixed to R = 10 8 , j m in = 10 and u\/c = 0.1. In the grayed 
region there is no stationnary state solution. 

We have reported in Fig. 6 the contour plot of a(l s , Ikin) and E c (l s ,lkin) for 
two values of e s (10 and 600 eV which are representative of the typical values 
of soft photons emitted by an accretion disk around a supermassive and stellar 
mass black hole respectively). The other parameters have been fixed to R = 10 8 , 
Imin = 10 and u\/c = 0.1. The contours of a and E c keep roughly the same 
shape but cover a different region of the parameter space for different parameter 
sets. In each figure, we can see that: 

• The spectral index does not strongly vary between the harder and the 
softer spectra (Act ~ 0.3). It reaches an asymptotic plateau for low l s and 
high l^in- In these conditions, both pair density and high energy cut-off 
are large so that the pair creation process is saturated. Concerning E c , 
following Eqs. (3) and (16), it is inversely proportional to l 2 s . 



12 



Petrucci et al. 



• The harder spectra are obtained for large values of l s or small values of 
Ikin- In both case, the pair efficiency decreases either due to a small value 
of E c or a low particle density, that is a low pair production optical depth. 
To keep the equilibrium between the pair creation effects and the hydro- 
dynamic of the flow, the system has to reach harder spectra to compensate 
this decrease of the pair creation rate. 

• For small l s , 7 C is very large (S> jth) an d its value becomes immaterial. 
Consequently, II and r, and thus a become independent of l s 

• For large lki n , the pair density must be also high to efficiently modify 
the hydrodynamical profile. The pair creation process is thus saturated 
meaning that the pair creation rate grows linearly with lki n , i-e. II, r and 
thus a become independent of lkin- 

• For to large l s or to small lkin, the system cannot reach sufficient hard 
states to keep in equilibrium and no high pair density stationary state can 
exist any more. The system can only be in the trivial pair free state (i.e. 
r = 7 and n = 0). 

We interpret the differences between the two plots of Fig. 6 as follows. The 
increase of e s favors the pair creation process. Thus, for given values of l s and 
lkin-, the spectral index is larger. The high energy cut-off E c increases mainly 
because of its dependence on e s (cf. Eq. (16)), 7 C keeping roughly constant (cf. 
Eq. (3)). 

5.3. Annihilation line 

The presence of pairs should give a signature as an annihilation feature at ~ 511 
keV. We will show that this feature is not expected to be strong in our model. 
Here, we have supposed the existence in the shock region of pre-accelerating pro- 
cesses bringing leptons to the sufficient energy (i.e. 7 > 7mm with 7mm 

of the 

order of a few, cf section 2.2) for resonant scattering off magnetic disturbances. 
Since the annihilation process occurs mainly at low energy, i.e. for particles with 
Lorentz factor 7 ~ 1 (cf. Coppi &: Blandford 1990), it occurs mainly far down- 
stream, where the pairs created in the shock can cool down. The annihilation 
line luminosity is thus at most equal to the pair rest mass luminosity. As shown 
in section 4, the pair luminosity II is itself limited and is necessarily smaller 
than ~20% of the X-ray/7-ray luminosity , i.e. II < H ma x — 0.2 (assuming 
that the total kinetic energy of the upstream flow is transformed in radiation). 
Besides, for II ~ 0.2, the compression factor is very small, of the order of unity, 
resulting in a very steep X-ray spectra. When hydrodynamics feedback is taken 
into account, the pair luminosity may be well below this theoretical limit of 
20%. An X-ray spectral photon index ~ 2 (as those generally seen in Seyfert 
galaxies) requires a compression ratio r ~ 3-4. Such values of r require values 
of n smaller than ~ 10% (cf. Fig 5). 

Assuming a steady state, pairs annihilate at the same rate as they are pro- 
duced. H/jmin gives then an upper limit of the annihilation line luminosity. 
We thus expect the luminosity of the annihilation radiation to be smaller than 
few percent of the total high energy radiation, which is quite compatible with 
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the non observation of strong annihilation lines in this class of Seyfert galaxies 
as shown by the best upper limit observed in Seyfert galaxies with the OSSE 
satellite (Johnson et al. 1997). 

5.4. Variability 

For some values of external parameters, as suggested by Fig. 6, no high pair 
density solution can exist in stationary states. So only the pair free solution 
(i,e, with II = 0) exists. The system is not expected to be variable with a 
constant set of parameters and variability can only occur with a variation of one 
of them. An interesting possibility would be to consider a possible feedback of 
the relativistic plasma to the soft compactness: in the reillumination models for 
instance (Collin, 1991; Henri & Petrucci, 1997), the soft photons are produced by 
the reprocessing of the primary X-ray emission. An increase of the pair plasma 
density will increase the X-ray illumination and thus the soft compactness. Fig. 
6 shows that in some cases, the change of l s make the system switch to pair free 
solution which will stop the reillumination and bring the system back to pair rich 
solutions. Limiting cycles could thus occur. We intend to further investigate 
this possible effect under astrophysically relevant conditions. 

6. Conclusion 

In the present paper, we have studied the effect of pair creation, via high energy 
photon-photon interaction, on a shock structure, where the high energy photons 
are produced via IC by the particles accelerated by the shock itself. The prob- 
lem is highly nonlinear since pairs can modify the shock profile through their 
pressure and, mutually, a change of the shock hydrodynamics can decrease or 
increase the pair production rate. 

We have shown that for a given size of the pair creation region, it exist 
a maximal value of the pair creation rate above which the shock cannot exist 
anymore. When the hydrodynamical feedbacks on the pair creation process are 
neglected, a pair power of at most 20% of the upstream kinetic power is sufficient 
to kill the shock. This constraint can fall to few percents in stationary states 
where pair creation and hydrodynamical effects balance. We thus do not expect 
the presence of strong annihilation lines. We also obtain spectral parameters in 
rough agreement with the observations. 

We suggest also a possible variability mechanism if the soft photon com- 
pactness depends itself on the pair density of the hot plasma, such as expected 
in reillumination models. 

In the model presented here, the cooling of particles is due to the IC process 
on external soft photons. However, particles may also cool on soft photon they 
produce by synchrotron process when spiraling around the magnetic field lines 
(the so-called synchrotron-self-compton process, SSC). In this case the cooling 
will also depend on the particle distribution function. We may expect that the 
addition of the SSC process would allow to obtain stationary states with harder 
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spectra than those we obtained here, since the additional synchrotron cooling 
would be compensate by a stronger acceleration, i.e. a larger compression ratio 
r. The detailed study of this problem is left to future work. 
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